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1 .  IimiODUCTION 


The  Prograa  Mane;;-!r  (PM)  for  BAT  (PM-BAT)  and  his  associated  contractors  (CAS, 
Incorporated,  ard  Dynetlcs,  Incorporated)  have  a  requirement  to  find  analytical 
distribution!-  to  adequately  represent  the  observed  wind  speed  distributions 
obtained  '  .  the  lower  S  km  of  the  atmosphere,  at  different  time  periods,  and  at 
different  geographical  regions  of  the  world.  These  distributions  are  required 
to  test  the  operational  capabilities  of  the  BAT  system.  The  U.S.  Army  Research 
Laboratory,  Battlefield  Environment  Directorate  (ARL-BED) ,  was  requested  to 
perform  the  task. 

Dramatic  changes  In  the  atmosphere  occur  each  day  on  the  average  as  the  sun  rises 
and  sets.  Because  diurnal  variations  are  most  prominent  in  the  atmosphere's 
lower  layer,  changes  in  boundary  layer  winds  have  been  documented  as  well  as  fer 
other  atmospheric  variables  (Holton,  1967).  Frequency  of  occurrence  of  low 
altitude  wind  speeds  has  important  applications  In  air  pollution  analysis  or  for 
nuclear-biological -chemical  (NBC)  predictions  for  military  applications.  The 
distributions  of  surface  wind  speeds  for  selected  time  periods,  geographical 
regions,  and  meteorological  conditions  have  been  reported  in  Avara  and  Miers 
(1992).  This  report  will  present  the  wind  speed  analysis  for  altitudes  above  the 
surface  -ip  to  S  km. 

The  upper-air  winds  are  obtained  by  Cracking  a  plbal  or  radiosonde  balloon  as  it 
ascends  Into  the  atmosphere.  The  relative  horizontal  movement  of  the  balloon  as 
it  traverses  successive  layers  of  the  atmosphere  can  be  determined  from  the 
positions  of  the  balloon  over  the  ground  at  successive  observing  times  as  It  Is 
tracked.  The  displacement  of  the  balloon  as  It  ascends  through  a  layer  is 
assumed  to  be  entirely  due  to  the  wind.  Furthermore,  Che  horizontal  velocity  of 
the  balloon  is  assumed  to  be  Che  horizontal  velocity  of  the  wind.  Wind  speed 
measurement  errors  typically  occur  through  balloon  tracking  eriors  and  srrors 
Induced  by  assuming  the  balloon  is  100  percent  responsive  to  the  wind  field. 
This  report  does  not  address  upper-air  wind  speed  measurement  errors.  The 
observed  wind  speeds  were  assumed  free  from  error. 

The  efforts  of  several  previous  Investigators  In  using  the  Ueibull,  gamma,  or 
log-normal  distributions  to  model  wind  speeds  are  reported  in  Avara  and  Miers 
(1992)  and  will  not  be  repeated  In  this  report.  All  three  distributions  were 
considered  in  Che  upper-air  wind  speed  analysis. 

2.  CLIMATOLOGY  DATA  BASE  AVAILABLE 

The  ARL-BED  maintains  a  climatology  data  base  containing  historical  upper-air 
weather  observations  from  numerous  locations  In  North  America,  Northern  and 
Central  Europe,  the  Mldeasc,  Southwest  Asia,  and  Korea.  These  data  were  obtained 
over  a  period  of  years  from  Che  U.S.  Air  Force  Environmental  Technical 
Applications  Center  (USAFETAC)  at  Scott  Air  Force  Base,  Illinois.  The  period  of 
record  for  the  data  from  each  station  is  generally  11  years  or  less  and  may  be 
anywhere  from  1966-76  to  1976-86,  depending  upon  the  year  the  data  were  obtained 
from  the  USAFETAC. 

The  regions  of  interest  were  the  European  Lowlands,  European  Highlands,  and  the 
Mldeasc  Desert  geographical  areas  (Miers  ec  al.,  1985a,  Miers  ec  al.,  1985b). 
The  upper- air  stations  located  in  each  of  these  regions  are  listed  In  table  1. 
ine  period  of  record  for  the  data  from  these  stations  Is  1973-83. 
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TABLE  1.  STATIONS  LOCATED  IN  THE  EDHOEEAN  LOWLANDS,  HIGHLANDS, 
AND  MIDEAS'I  DESE&T  CECCXAEHICAL  REGIONS 


Station 

Number 

Station 

Name 

Latitude 

(DegiMln) 

Longitude 

(Deg;Mln) 

Elevation 

(Meters) 

EUROPEAN  LOWLANDS 

100350 

Schleswig,  FRG 

+54:32 

-009:33 

+0048 

100460 

Klel/Holtenau,  FRG 

+54:23 

-010:09 

+0031 

102380 

Porgin/Hohne .  FRG 

+32:49 

-009:56 

+0068 

103040 

Meppen,  FRG 

+52:43 

-007:19 

+0019 

103380 

Hannover,  "RG 

+52:28 

-009:42 

+0055 

103840 

Borlin/Templehof ,  FRG 

+52:28 

-013:24 

+0050 

104100 

Essen,  FRG 

+51 ;24 

-006:58 

+0161 

EPROrgAP-HICHURP? 

095480 

Melnlngen,  DD 

+50:34 

-010:23 

+045C 

106180 

IDAR/Obersteln,  FRG 

+49:42 

-007:20 

+0376 

106870 

Graftenwohr,  FRG 

+49:42 

-011:57 

+0414 

107710 

Garmersdorf,  FRG 

+49:26 

•011:54 

+0419 

HIP.EAS.T  aSSgfxT 

400300 

Damascus  New  Inti,  SY 

+33:25 

-036:31 

+0610 

402650 

Mafraq,  TJ 

+32:22 

-036:16 

+0686 

406500 

Baghdad  Int  (ORBW) .  IQ 

+33:14 

-044:14 

+0034 

407540 

Tehran/Mehrabad,  IR 

+35:41 

-051:19 

+  1204 

For  the  two  European  regions,  the  winter  and  annual  wind  speed  distributions  were 
obtained  by  grouping  the  observations  fron  all  hours  of  the  day  together.  Winter 
for  these  regions  is  defined  to  be  November,  December,  January,  and  February. 
For  the  Hldeast  Desert  region,  the  summer  and  annual  wlna  speed  distributions 
were  obtained  by  grouping  the  observations  taken  at  nighc  into  one  group  and  the 
observations  taken  during  the  daytime  into  another  group.  Summer  for  the  Hldeast 
Desert  region  is  defined  to  be  June,  July,  and  Augjst. 

The  upper-air  observations  tend  to  be  taken  at  0000  and  1200  CHT  with  possibly 
a  few  observations  taken  at  other  times  of  day  on  a  nrnroutlne  basis.  Therefore, 
performing  an  analysis  of  upper-air  data  over  more  than  two  periods  of  the  day 
is  not  practical  for  any  geographical  region  or  season. 

The  historical  upper-air  wind  speed  observations  were  linearly  interpolated  to 
give  wind  speeds  every  500  m  from  500  m  altitude  up  to  5000  m.  All  altitudes 
referred  to  in  this  report  are  relative  to  mean  sea  level  (MSL1 .  The  wind  speeds 
wt  these  10  altitude  levels  were  used  In  the  analysis. 
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3.  ANALYTICAL  DISTRIBUTIONS 


Three  analytical  distributions  were  selected  for  potential  use  in  representing 
the  observed  wind  speed  distributions.  They  were  chosen  based  upon  the  work  of 
previous  researchers  referred  to  In  Avara  and  Mlers  (1992).  These  three 
candidate  distributions  were  the  ganuaa,  log-normal,  and  Weibull,  Each  of  these 
distributions  performed  reasonably  well  at  representing  the  observed  wind  speed 
distributions . 

The  gamma  probability  distribution  is  defined  by 


f  P) 


(-i)<*-‘'exp[-(-^)3 


where  r(a)  is  the  gamma  function  with  argument  a,  x  >  0.  o  >  0.  and  ^  >  0. 
The  iog-normal  probability  distribution  is  defined  by 


f(x;\k.a)  • 


exp(-4(-i5i£L:li)^] 
_ 2 _ _ _ 


where  n  is  the  expected  value  and  a  is  the  standard  deviation  of  ln(x)  and  x  >  0, 
The  Weibull  probability  distribution  is  defined  by 


=  (|)  (-^)-’exp(-(-^)']  , 


where  x  >  0,  a  >  0,  and  ^  >  0. 

In  this  report  p  and  ^  are  referred  to  as  scale  parameters  because  changing  their 
values  either  tends  to  move  the  center  or  skew  the  distribution  toward  higher  or 
lower  values,  a  and  a  are  referred  to  as  shape  parameters  because  changing  their 
values  tends  to  change  the  shape  of  the  distribution  to  either  a  more  peaked  or 
a  flatter  form. 

Three  procedures  for  computing  values  for  a  and  fi  ox  n  and  a  were  considered  for 
fitting  che  analytical  distributions  to  the  observed  wind  speed  distributions 
(histograms).  These  procedures  are  the  method  of  moments,  maximum  likelihood, 
and  nonlinear  least  squares  fit  to  the  observed  cumulative  probability  distribu¬ 
tion.  For  each  observed  wind  speed  distribution  the  best  of  the  three  analytical 
distributions  and  the  associated  parameter  values  were  determined.  Once  the  name 
of  che  best  analytical  distribution  and  che  distribution  parameter  values  are 
known,  the  estimation  of  the  Call  probability  of  che  wind  speed  is  possible. 
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The  “goodness-of -  fit*  of  the  analytical  distributions  was  compared  by  computing 
both  Che  mean  absolute  difference  and  the  rooc-mean-square  (rms)  difference 
between  the  observed  and  each  of  the  three  analytical  cumulative  distributions 
After  thoroughly  comparing  the  analytical  distributions  by  using  bt'fh  criteria, 
the  author  decided  to  use  the  minimum  rms  difference  to  make  the  final  sel.  -ion. 

When  this  criterion  is  used,  the  method  of  moments  generally  yields  parameter 
estimates  that  do  not  perform  as  well  as  Che  other  two  procedures.  The  maximum 
likelihood  procedure  tends  to  yield  better  estimates  of  the  pe-,.an!ecers  than  the 
method  of  moments.  The  nonlinear  least  squares  fit  giver  the  best  parameter 
estimates,  as  exp^'cted  by  the  selection  criterion. 

3.1  Method  of  Momenta 

The  expected  value  of  x  for  the  gamma  distribution  is  given  by 

EtxJ  -  ttP,  (*) 


and  the  variance  of  x  is  given  by 

VAR[.x]  »  ap^.  ^5) 


These  two  equations  can  be  used  Co  compute  values  of  a  and  ^  for  the  gamma 
distribution  using  observed  data  to  estimate  E[x]  by  the  mean  value  of  the 
observed  data  and  VAR(xj  by  Che  sample  variance  of  the  observed  data,  giving 

a  ,  VAHlxl  (6) 


and 

a  =  ,  (7) 

P 


This  procedure  for  obtaining  values  for  the  distribution  parameters  is  called  the 
method  of  moments. 

Using  the  method  of  moments  Co  compute  values  for  n  and  a  in  the  log-norraal 
distribution  gives 

p  -  CllnCx) ]  (8) 
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and 


o  ■  VVAR I in (xj  j , 


(9) 


where  E[ln<x)]  and  VAR(ln(x)J  are  eaclaaced  froa  the  daca  by  using  the  aean  value 
and  sample  variance  o£  the  natural  logaricha  of  the  wind  speed  data. 

The  advantage  of  using  E[ln(x)l  and  VAR(ln(x)]  rather  than  E(x]  and  VAR(xl  to 
estimate  n  and  c  can  be  seen  when  one  considers  the  relationship  between  the 
distribution  parameters  and  E[x]  and  VAR[x]  given  by 

ECx]  ■  exp[|»*-jO*3 


and 


VASlx]  ■  exp(2p*o*)  (exp{o*)  -1]  . 


(11) 


The  expected  value  of  x  for  the  Welbull  distribution  is  given  by 

Fix]  »  |r(l). 

8 


(12) 


and  the  variance  of  x  la  given  by 


varix]  -  i!(2r(l)  -  Ain-i)]*]. 

8  a  8  a 


(13) 


It  is  rather  difficult  to  use  the  method  of  moments  to  compute  values  of  a  and 
P  for  the  Welbull  distribution  using  estimates  of  Elx]  and  VAR[x]  obtained  from 
observed  data.  When  the  computational  time  required  to  iteratively  solve  this 
system  of  equations  and  the  general  inferior  quality  of  the  parameter  estimates 
were  considered,  the  author  did  not  use  this  procedure  to  compute  the  Welbull 
distribution  parameters. 

3.2  Maximum  Likelihood 

The  maximum  likelihood  procedure  Involves  finding  the  parameter  values  that 
maximize  the  likelihood  function  of  the  distribution.  For  the  gaiama  distribution 
these  parameter  values  are  obtained  by  iteratively  solving  the  following  system 
of  equations: 
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ln<P>  ♦  ¥(<»)  -  S[ln(x)] , 


(14) 


af  -  £ix] , 


(15) 


where  t(a)  is  the  dlgamms  function  defined  by  the  equation 

-  dfln(r(aU) 


(1€) 


with  r(a)  being  the  ganoa  function  with  arguiaent  a  (Abranowltz  and  Stegun,  1965) . 

The  oaximua  likelihood  estimates  of  the  log-normal  distribution  parameters  are 
the  same  as  chose  given  in  equations  (8)  and  (9)  by  Che  method  of  moments. 

The  maximum  likelihood  estimates  of  the  Weibull  distribution  parameters  are 
obtained  by  iteratively  solving  the  following  system  of  equations: 


g(x‘ln(x)l  _  i 
Blx*]  « 


Sfln(x) ) , 


(17) 


1 

P  -  K{X»)  •  . 


(18) 


3.3  Nonlinear  Least  Squares  Fit 

The  least  squares  fit  was  performed  on  observed  cumulative  distribution  values 
(Pj,  1  -  1,  2,  ...  ,  NJ  corresponding  to  discrete  values  of  the  wind  speed,  X, 
denoted  by  (X^,  1  -  1,  2,  ...  ,  N)  where  H  is  the  number  of  wind  speed  intervals. 
Each  wind  speed  interval  had  a  width  of  0.5  m/s.  N  was  determined  by  choosing 
the  first  interval  of  the  observed  wind  speed  frequency  of  occurrence  histogram 
such  that  all  intervals  greater  than  the  Nth  contained  no  occurrences  of  wind 
speed. 

The  three  analytical  distributions  are  clearly  nonlinear  with  respect  to  Che 
distribution  parameters  a,  n,  and  a.  To  determine  Che  least  squares  estimates 
of  these  parameter  values,  a  procedure  was  used  Co  linearize  the  distributions 
(at  least  in  the  neighborhood  of  any  specific  values  of  the  parameters). 
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Consider  the  integral  of  the  probability  distribution  given  by 


f’<X;a.P)  •  ff(x;a.P)dx. 


(19) 


By  letting 


a  •  Sg  *  Aa 


(20) 


and 


P  -  P,  ♦  89 


(21) 


equation  (19)  becooes 


FiXiu^  *  Aa.p,  ♦  AP)  -  ♦  Aa.P,  ♦  Ap)dx. 


(22) 


To  a  first'order  approxination  using  snail  values  for  Sa  and  30 


FlX;a,fi}  -f’(^;a,,Pg)  ♦ 


af'(Jf;a„pg)  ^  af(X;a,.pg) 

- 35 - *  - 53 - 


AP, 


(23) 


where 


af'(X;aa.  Pg)  _  r  df(x!a,  8)  i  ^ 

35  J  da  .J., 


(24) 


and 


dF(X;ag,Pg)  _  fdf(x;a,B) 

— 33 - J— -iJ.. 


dx. 


(25) 


Equation  (23)  is  linear  in  da  and  30. 
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Let  Q  represent  the  sua  of  the  squerea  of  the  differences  between  the  enelytlcal, 
F(Xi;  a,  fi) ,  and  observed  cuaulatlve.  Pj.  distributions  given  by 


(26) 


Replacing  F(Xj;  a,  fi)  by  the  expression  in  equation  (23)  gives 


(27) 


Taking  the  partial  derivatives  of  Q  with  respect  to  6a  and  6fi,  and  setting  each 
derivative  equal  to  zero  to  find  a  ninioua  of  Q  gives 


and 


^  ,iw iw ^  ^  ,  «,> 


where  the  terns  involving  F(](Xi)  and  its  partial  derivatives  denote  the  function 
F(x;  a,  0)  and  its  partial  derivatives  evaluated  at  X^,  oo,  and  0o.  These  two 
eqxiatlons  are  linear  In  6a  and  60. 

The  nonlinear  least  squares  fit  procedure  begins  by  getting  Initial  values  for 
og  and  0g;  solving  equations  (28)  and  (29)  for  da  and  60  using  the  definitions 
in  equations  (19),  (2A),  and  (25);  and  getting  updated  values  of  a  and  0 
according  to  the  iterative  equations 


«K  •  a,^.,  ♦  8a  (30) 

and 

P*  -  P.-,  ♦  «P  (31) 


for  k  -  1,  2,  ...  until  6a  and  60  becone  very  snail  and  convergence  is  obtained. 
The  initial  values  of  the  paraaeters  nay  be  selected  by  any  nethod.  Reasonable 
initial  values  are  obtained  by  using  the  nethod  of  nonents  to  get  Oq  and  0g. 
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4.  tf^ND  SPEED  DISTRIBUTIONS 


All  three  of  the  analytical  distributions  considered  in  this  analysis  have  zero 
probability  of  occurrence  at  x  •  0.  The  observed  wind  speed  at  any  altitude, 
however,  may  be  cairn  or  zero.  To  account  for  this  phenomenon,  the  wind  speed 
cumulative  distribution  is  defined  by  the  equation 


W(.Jf;Po,«,P)  -p,8<0)  ♦  (1  -  Po)f(X;a,p)  , 


(32) 


where  X  is  the  wind  speed,  d(0>  is  the  Dirac  delta  function,  pg  is  the 
probability  of  calm  wind  speeds,  and  F(X;  a,  /3)  is  the  appropriate  analytical 
cumulative  distribution  function. 

The  probabilities  of  calm  wind  speeds  for  each  region,  time  period,  and  height 
level  combination  are  listed  in  table  2.  Similarly,  the  numbers  of  noncalm  wind 
speed  observations  are  listed  in  table  3. 

The  observed  probability  distribution  was  obtained  by  forming  histograms  of 
frequency  of  occurrence  of  observed  wind  speed  with  all  histogram  Interval  widths 
being  0.5  m/s.  The  histogram  frequencies  of  occurrence  were  converted  to  proba¬ 
bilities  and  assigned  to  the  wind  speeds  at  the  midpoints  of  the  intervals.  The 
corresponding  analytical  probability  distribution  was  obtained  by  evaluating  the 
analytical  cumulative  probability  distribution  at  the  wind  speeds  corresponding 
to  Che  beginning  and  end  of  each  histogram  interval  and  using  the  difference 
between  the  two  values  assigned  to  the  wind  speed  at  Che  midpoint  of  Che 
interval . 

The  mean  absolute  deviation  between  the  observed  and  analytical  cumulative 
distributions  for  each  region,  time  period,  and  height  level  combination  is 
listed  in  cable  4.  Similarly,  the  rms  deviation  between  the  cumulative  distri¬ 
butions  is  listed  in  table  5.  Generally,  both  the  mean  absolute  deviations  and 
the  rms  deviations  are  small,  indicating  the  analytical  distributions  adequately 
represent  the  observed  wind  speed  distributions.  The  names  of  the  corresponding 
best  distribution,  Che  shape  parameters,  and  the  scale  parameters  are  listed  in 
Cables  6,  7,  and  8,  respectively. 
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TABLE  2 .  PROBABILITIES  OF  CALX  WIND  SFEEOS  FOR  ALL  WIND  SPEED  DISTRIBUTION  CASES 


Height 


_ Season 

Tine 

Period 

Level 

fkia^ 

Probability 

fXI 

European  lowlands 

Winter 

All  hours 

0.5 

.09317 

1.0 

.04192 

1.5 

.02474 

2.0 

.01651 

2.5 

.02064 

3.0 

.04132 

3.5 

.02481 

4.0 

.01655 

4.5 

.01672 

5.0 

.00836 

European  lowlands 

Annual 

All  hours 

0.5 

.12822 

1.0 

.05535 

1.5 

.04233 

2.0 

.03143 

2.5 

.03145 

3.0 

.05062 

3.5 

.02465 

4.0 

.02192 

4.5 

.01665 

5.0 

.OHIO 

European  highlands 

Winter 

All  hours 

0.5 

.50052 

1.0 

.05432 

1.5 

.02647 

2.0 

.00000 

2.5 

.00000 

3.0 

. 00000 

3.5 

.00000 

4.0 

. 00000 

4.5 

.00000 

5.0 

,00000 

European  highlands 

Annual 

All  hours 

0.5 

.62929 

1.0 

.13030 

1.5 

.06565 

2.0 

.02195 

2.5 

.02642 

3.0 

.03093 

3.5 

. 00444 

4.0 

.00445 

4.5 

.00000 

5.0 

.00450 
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Mideast  desert 

Sutamer 

Night 

0.5 

1.59272 

1.0 

.74587 

1.5 

.33656 

2.0 

.12573 

2.5 

.12658 

3.0 

.08460 

3.5 

. 00000 

4.0 

.00000 

4.5 

. 04456 

5.0 

.04486 

Hldeast  desert 

Sunmer 

Day 

0.5 

.93555 

1.0 

.35858 

1.5 

.06165 

2.0 

.06094 

2.5 

.09146 

3.0 

.15263 

3.5 

.06139 

4.0 

.03095 

4.5 

.00000 

5.0 

.00000 

Mideast  desert 

Annual 

Night 

0.5 

3.26118 

1.0 

1.15770 

1.5 

.56747 

2.0 

.08887 

2.5 

.05568 

3.0 

.04468 

3.5 

.00000 

4.0 

.00000 

4.5 

.01177 

5.0 

.01188 

Mideast  desert 

Annual 

Day 

0.5 

2.90059 

1.0 

.95894 

1.5 

.49075 

2.0 

.12841 

2.5 

.07709 

3.0 

.08570 

3.5 

.01725 

4.0 

.01735 

4.5 

.00000 

5.0 

.00000 
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TABLE  3.  THE  RUMIEE  OF  HONCAXM  VIND  SPEED  OBSERVATIONS  FOR  ALL  WIND  SPEED 
DISTRIBUTION  CASES 


Revlon 

Season 

Time 

Period 

Height 

Level 

Noncalm 

Count 

European  lowlands 

Winter 

All  hours 

0.5 

23592 

1.0 

23845 

1.5 

24251 

2.0 

24230 

2.5 

24215 

3.0 

24194 

3.5 

24182 

4.0 

24170 

4,5 

23918 

5.0 

23913 

European  lowlands 

Annual 

All  hours 

0.5 

71661 

1.0 

72230 

1.5 

73202 

2.0 

73157 

2.5 

73106 

3.0 

72056 

3.5 

73017 

4.0 

72983 

4.5 

72066 

5.0 

72055 

European  highlands 

Winter 

All  hours 

0.5 

6759 

1.0 

7360 

1.5 

7555 

2.0 

7547 

2.5 

7530 

3.0 

7508 

3.5 

7503 

4.0 

7478 

4.5 

7438 

5.0 

7395 

European  highlands 

Aimual 

All  hours 

0.5 

20686 

1.0 

22227 

1.5 

22832 

2.0 

22773 

2.5 

22701 

3.0 

22623 

3.5 

22540 

4.0 

22471 

4.5 

22343 

5.0 

22220 

Kldeast  desert 

Suoner 

Night 

0.5 

865 

1.0 

1863 

1.5 

2369 

2.0 

2383 

2.5 

2367 

3.0 

2362 

3.5 

2304 

4.0 

2268 

4.5 

2243 

5.0 

2228 

Mideesc  desert 

Suoaer 

Day 

0.5 

953 

1.0 

2223 

1.5 

3242 

2.0 

3280 

2.5 

3277 

3.0 

3271 

3.5 

3256 

4.0 

3230 

4.5 

3208 

5.0 

3192 

HldeasC  desert 

Annual 

Night 

0.5 

3352 

1.0 

7001 

1.5 

8761 

2.0 

8994 

2.5 

8975 

3.0 

8949 

3.5 

8734 

4.0 

8614 

4.5 

8496 

5.0 

8416 

Mideast  desert 

Annual  Day 

0.5 

3448 

1.0 

8036 

1.5 

11355 

2.0 

11666 

2.5 

11666 

3.0 

11658 

3.5 

11590 

4.0 

11524 

4.5 

11454 

5.0 

11397 
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TABU  4.  THE  HEAM  ABSOLUTE  DEVIATION  BETWEEN  THE  OBSEEVSO  AND  ANALmCAL 
CUMULATIVE  DISTRIBUTIONS 


Re£ion 

Season 

Tiaa 

Period 

Height 

Level 

(kn) 

Mean  Absolute 
Deviation 

(X) 

European  lowlands 

Winter 

All  hours 

0.5 

.14902 

1.0 

.28980 

1.5 

.32437 

2.0 

.26155 

2.5 

.26858 

3.0 

.24476 

3.5 

.23916 

4.0 

.16868 

4.5 

.11883 

5.0 

.09139 

European  lowlands 

Annual 

All  hours 

0.5 

.20985 

1.0 

.29884 

1.5 

.29448 

2.0 

.24106 

2.5 

.26972 

3.0 

.29805 

3.5 

.30699 

4.0 

.27436 

4.3 

.25471 

5.0 

.23513 

European  highlands 

Winter 

All  hours 

0.5 

.38035 

1.0 

.30137 

1.5 

.30301 

2.0 

.40481 

2.5 

.28888 

3.0 

.19688 

3.5 

.18850 

4.0 

.25313 

4.5 

.24869 

5.0 

.20566 

European  highlands 

Annual 

All  hours 

0.5 

.42646 

1.0 

.26195 

1.5 

.21188 

2.0 

.12471 

2.5 

.22099 

3.0 

.15987 

3.5 

.16872 

4.0 

.20203 

4.5 

.21563 

5.0 

.16645 

Mlduaat  desert  Suoaer  NlghC  0.5  1.07418 

1.0  .38178 

1.5  .65171 

2.0  .29377 

2.5  .39364 

3.0  .57006 

3.5  .44411 

4.0  .35257 

4.5  .28585 

5.0  .29915 

Mideast  desert  Sumaer  Day  0.5  .30531 

1.0  .51226 

1.5  .45835 

2.0  .30314 

2.5  .34658 

3.0  .22126 

3.5  .49947 

4.0  .35426 

4.5  .31233 

5.0  .30977 


Mideast  desert  Annual  Night  0.5  .31784 

1.0  .23407 

1.5  .56314 

2.0  .25396 

2.5  .14807 

3.0  .36652 

3.5  .12457 

4.0  .20725 

4.5  .12137 

5.0  .14451 


Mideast  desert  Annual  Day  0.5  .43033 

1.0  .25526 

1.5  .26052 

2.0  .24431 

2.5  .12826 

3.0  .23874 

3.5  .09874 

4.0  .25631 

4.5  .32095 

5.0  .25691 
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TABLE  S.  TUZ  lOOT-XBAN-SQUAU  DEVUTIOS  BETVEEB  THE  OBSERVED  AND  ANALYTICAL 
CDHDUTIVE  OISTKIBDTIONS 


Ration 

Season 

TIm 

Period 

Height 

Level 

fka) 

Root  .Kean  Squj 
Deviation 
ai 

Europatn  lowlands 

Winter 

All  Hours 

0.5 

.29728 

1.0 

.44189 

1.5 

.50877 

2.0 

.42723 

2.5 

.43830 

3.0 

.37937 

3.5 

.35263 

4.0 

.26514 

4.5 

.21340 

5.0 

.18355 

European  lowlands 

Annual 

All  Hours 

0.5 

.30101 

1.0 

.46328 

1.5 

.4499^ 

2.0 

.36094 

2.5 

. 34462 

3.0 

.38526 

3.5 

.40791 

4.0 

.42809 

4.5 

.40682 

5.0 

.40449 

European  highlands 

Winter 

All  Hours 

0.5 

.61938 

1.0 

.42085 

1.5 

.46469 

2.0 

.58212 

2.5 

.45906 

3.0 

.30152 

3  5 

.29192 

4.0 

.35293 

4.5 

.34589 

5.0 

.29311 

European  highlands 

Annual 

All  Hours 

0.5 

.79256 

1.0 

.32503 

1.5 

.31076 

2.0 

.22408 

2.5 

.37384 

3.0 

.34483 

3.5 

.36562 

4.0 

.39665 

4.5 

.41282 

5,0 

.33166 
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Mideast  desert 

Suaaer 

Night 

0.5 

1.24607 

1.0 

.55307 

1.5 

.80162 

2.0 

.42284 

2.5 

.53262 

3.0 

.73705 

3.5 

.64822 

4.0 

.48737 

4.5 

.43040 

5.0 

.48113 

Mideast  desert 

Sumei 

Day 

0.5 

.42728 

1.0 

.83192 

1.5 

.67613 

2.0 

.48374 

2.5 

.50753 

3.0 

.30667 

3.5 

.69918 

4.0 

.51209 

4.5 

.43699 

5.0 

.48574 

Mideast  desert 

Annual 

Night 

0.5 

.47518 

1.0 

.44673 

1.5 

.85175 

2.0 

.36001 

2.5 

.21746 

3.0 

.58937 

3.5 

.26273 

4.0 

.35700 

4.5 

.27835 

5  0 

.34835 

Mideast  desert 

Annual 

Day 

0.5 

.60093 

1.0 

.46991 

1.5 

.48737 

2.0 

.44501 

2.5 

.21835 

3.0 

.38982 

3.5 

.15114 

4.0 

.39980 

4.5 

.48940 

5.0 

.39138 
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TABLE  «.  THE  NAXES  OF  THE  BEST  ANALTTICAL  DISTRIBUTIONS  TO  REFRESENT  THE 
OBSERVED  VIND  SPEED  DISTRIBUTIONS 


European  lowlands  Winter  All  hours  0.5  Uelbull 

l.O  Welbull 

1.5  Welbull 

2.0  Welbull 

2.5  Welbull 

3.0  Welbull 

3.5  Welbull 

4.0  Welbull 

4.5  Welbull 

5.0  Welbull 


European  lowlands  Annual  All  hours  0 . 5  Gamoa 

1.0  Garama 

1 . 5  Canffla 

2 . 0  Ganna 

2 . 5  Ganna 

3 . 0  Ganna 

3 . 5  Ganna 

4.0  Welbull 

4.5  Welbull 

5.0  Welbull 


European  highlands  Winter  All  hours  0.5  Ganna 

1.0  Welbull 

1.5  Welbull 

2.0  Welbull 

2.5  Welbull 

3.0  Welbull 

3.5  Welbull 

4.0  Welbull 

4.5  Welbull 

5.0  Uelbull 


European  highlands  Annual  All  hours  0.5  Uelbull 

1 . 0  Ganna 

1 . 5  Ganna 

2.0  Ganna 

2 . 5  Ganna 

3.0  Welbull 

3.5  Welbull 

4.0  Welbull 

4.5  Welbull 

5.0  Welbull 
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Mideast  desert 

Suamer 

Night 

0.5 

Gama 

1.0 

Gama 

1.5 

CaoBia 

2.0 

Gama 

2.5 

Gama 

3.0 

Log'Nomal 

3,5 

Gama 

4.0 

Gama 

4.5 

Gama 

5.0 

Gama 

Mideast  desert 

Sumer 

Day 

0.5 

Weibull 

1.0 

tfelbull 

1.5 

Gama 

2.0 

Gama 

2.5 

Gama 

3.0 

Gama 

3.5 

Gamaa 

4.0 

GaoBBa 

4.5 

Gama 

5.0 

Gama 

Mideast  desert 

Annual 

Night 

0.5 

Gama 

1.0 

Gama 

1.5 

Weibull 

2.0 

Gama 

2.5 

Gama 

3.0 

Gama 

3.5 

Gama 

4.0 

Gajraa 

4.5 

Gama 

5.0 

Gama 

Mideast  desert 

Annual 

Day 

0.5 

Weibull 

1.0 

Ganna 

1.5 

Ganna 

2.0 

Gama 

2.5 

Ganna 

3.0 

Gama 

3.5 

Gama 

4.0 

Gama 

4.5 

Gama 

5.0 

Ganna 
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TABUt  7 .  THE  SHAPE  PARAMETER  CORRESPONDIMO  TO  THE  BEST  AHALYTICAL  DISTRIEDTION 
FOR  ALL  VINO  SPEED  DISTRIBUTION  CASES 


Height 


Time 

Period 

Level 

(to) 

Shape 

Paramete; 

European  lowlanda 

Winter 

All  hours 

0.5 

1.94001 

1.0 

1.87258 

1.5 

1.87244 

2.0 

1.90308 

2.5 

1.93556 

3.0 

1.92341 

3.5 

1.95354 

4.0 

1.93692 

4.5 

1.93676 

5.0 

1.92517 

European  lowlands 

Annual 

All  hours 

0.5 

2.99105 

1.0 

2.93947 

1.5 

2.94049 

2.0 

3.05194 

2.5 

3.03083 

3.0 

2.88379 

3.5 

3.00621 

4.0 

1.83622 

4.5 

1.82808 

5.0 

1.81675 

European  highlands 

Winter 

All  hours 

0.5 

2.04370 

1.0 

1.74782 

1.5 

1.72693 

2.0 

1.74782 

2.5 

1.78564 

3.0 

1.79718 

3.5 

1.83906 

4.0 

1.85042 

4.5 

1.86825 

5.0 

1.86339 

European  highlands 

Annual 

All  hours 

0.5 

1.58353 

1.0 

2.59101 

1.5 

2.49970 

2.0 

2.58219 

2.5 

2.72753 

3.0 

1.73934 

3.5 

1.76690 

4.0 

1.77797 

4.5 

1.78218 

5.0 

1.77553 
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Hldeast  desert 

SuBBier 

Might 

0.5 

3.54770 

1.0 

2.97721 

1.5 

3.06661 

2.0 

3.76824 

2.5 

3.98637 

3.0 

0.54431 

3.5 

3.44719 

4.0 

4.01631 

4.5 

3.98119 

5.0 

3.98657 

HideasC  desert 

Suaner 

Day 

0.5 

2.74135 

1.0 

1.95323 

1.5 

3.42313 

2.0 

3.93455 

2.5 

3.92779 

3.0 

3.34170 

3.5 

3.89361 

4.0 

4.04577 

4.5 

4.13176 

5.0 

3.98775 

Mideast  desert 

Annual 

Might 

0.5 

3.00832 

1.0 

2.48724 

1.5 

1.80496 

2.0 

3.28748 

2.5 

3.51753 

3.0 

3.14780 

3.5 

3.45259 

4.0 

3.65231 

4.5 

3.81167 

5.0 

3.75729 

Mideast  desert 

Annual 

Day 

0.5 

2.21271 

1.0 

2.95545 

1.5 

2.81773 

2.0 

3.40679 

2.5 

3.49887 

3.0 

3.12657 

3.5 

3.47248 

4.0 

3.68069 

4.5 

3.73394 

5.0 

3.61816 
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TABLE  8.  THE  SCALE  PAXAMETEE  CORRESFOKCIHC  TO  THE  BEST  AMALTTICAL  DISTRIBUTION 
FOR  ALL  WIND  SPEED  DISTRIBUTION  CASES 


Helghc 


Tiow 

Level 

Scale 

European  lowlands 

Winter 

All  hours 

0.5 

10.52280 

1.0 

13.23674 

1.5 

13.76583 

2.0 

14.18878 

2.5 

14.81256 

3.0 

15.54027 

3.5 

16.56613 

4.0 

17.62310 

4.5 

18.76522 

5.0 

19.90936 

European  lowlands 

Annual 

All  hours 

0.5 

2.63915 

1.0 

3.25972 

1.5 

3.42868 

2.0 

3.56763 

2.5 

3.74912 

3.0 

3.95533 

3.5 

4.21879 

4.0 

14.91398 

4.5 

15.87218 

5.0 

16.84511 

European  Highlands 

Winter 

All  Hours 

0.5 

2.29426 

1.0 

10.29712 

1.5 

12.23363 

2.0 

13.03796 

2.5 

13.72103 

3.0 

14.52975 

3.5 

15.59311 

4.0 

16.76794 

4.5 

17.92975 

5.0 

19.10730 

European  highlands 

Annual 

All  hours 

0.5 

4.55167 

1.0 

3.07453 

1.5 

3.63359 

2.0 

3.92955 

2.5 

4.18699 

3.0 

12.35809 

3.5 

13.28000 

4.0 

14.31750 

4.5 

15.33360 

5.0 

16.34973 
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Midflsst  desert 

SuoBser 

Might 

0.5 

2.33842 

1.0 

2.19361 

1.5 

2.24858 

2.0 

1.72230 

2.5 

1.78385 

3.0 

1.89340 

3.5 

2.33097 

4.0 

2.16897 

4.5 

2.34215 

5.0 

2.52840 

Mldeesc  desert 

Sumer 

Day 

0.5 

9.18078 

1.0 

7.42801 

1.5 

1.82355 

2.0 

1.62747 

2.5 

1.70796 

3.0 

2.07340 

3.5 

2.01106 

4.0 

2.17162 

4.5 

2.35531 

5.0 

2.51934 

Mideast  desert 

Annual 

Night 

0.5 

2.26933 

1  0 

2.30821 

1.5 

7.11873 

2.0 

2.03827 

2.5 

2.29445 

3.0 

3.01649 

3.5 

2.87908 

4.0 

3.21601 

4.5 

3.08457 

5.0 

3.37354 

Hldeast  desert 

Annual 

Day 

0.5 

7.53769 

1.0 

2,01991 

1.5 

2 .  ■'  1044 

2.0 

1.99963 

2.5 

2.21907 

3.0 

2.91417 

3.5 

2.82217 

4.0 

3.16577 

4.5 

3.49970 

5.0 

3.82600 
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5 .  ANALYSIS  RISULTS 


5 . 1  European  Lowlands 

The  observed  noncalm  wind  speed  and  analytical  distributions  along  with  the 
observed  and  analytical  cumulative  distributions  for  the  European  Lowlands  are 
depicted  in  figures  1  through  20.  The  figures  correspond  to  winter  and  annual 
time  periods  for  the  10  height  levels.  The  upper  part  of  each  figure  depicts 
the  observed  and  analytical  distributions  while  the  lower  part  of  the  figure 
depicts  the  observed  and  analytical  cumulative  distributions.  In  each  figure  the 
differences  between  the  observed  and  analytical  values  are  also  plotted.  The 
observed  distributions  are  indicated  by  a  dashed  line,  the  analytical  distri¬ 
butions  by  a  solid  line,  and  the  difference  between  them  by  a  dotted  line. 

The  figures  depict  only  the  noncalm  part  of  the  wind  speed  distributions.  The 
number  of  noncalm  wind  observations  used  to  fit  the  analytical  distributions  to 
the  observed  wind  speed  distributions  Is  displayed  in  each  figure.  In  addition, 
the  best  analytical  probability  distribution  and  the  corresponding  scale  and 
shape  parameters  are  listed  In  each  figure. 

Of  the  20  different  height  level  and  season  combinations  there  were  13  cases 
where  the  Welbull  distribution  was  the  best  and  7  where  the  gamma  was  the  best 
(reference  table  6).  The  Welbull  distribution  was  the  best  for  all  winter  wind 
speed  distributions  and  for  annual  wind  speed  distributions  for  height  levels 
above  3.5  km.  The  gamma  distribution  was  the  best  for  annual  wind  speed  distri¬ 
butions  for  height  levels  up  to  3.5  km.  Overall,  the  analytical  distributions 
perform  well  at  describing  the  observed  wind  speed  distributions  for  the  European 
lowlands . 

5 . 2  European  Highlands 

Figures  21  through  40  depict  tha  same  type  of  wind  speed  information  for  the 
European  highlands.  Of  Che  20  different  height  level  and  season  combinations 
there  were  15  cases  where  the  Welbull  distribution  was  the  best  and  5  where  the 
gamma  was  the  best  (reference  table  6).  The  Welbull  distribution  was  best  for 
winter  wind  speed  distributions  for  height  levels  above  0.5  km  and  for  all  annual 
wind  speed  distributions  except  those  with  height  levels  greater  than  0.5  km  and 
less  than  3.0  km.  The  gamma  distribution  was  the  best  for  the  winter  wind  speed 
distribution  at  height  level  0.5  km  and  Che  annual  wind  speed  distributions  for 
height  levels  greater  than  0.5  km  and  less  than  3.0  km.  Overall,  the  analytical 
distributions  perform  well  at  describing  the  observed  wind  speed  distributions 
for  the  European  highlands. 

5.3  Mldeast  Desert 

Figures  41  through  80  depict  the  same  type  of  wind  speed  information  for  the 
Mldeast  desert  region  except  that  the  distributions  are  for  summer  and  annual  for 
nighttime  and  daytime  periods.  Of  the  40  different  height  level  and  time  period 
combinations  there  were  35  cases  where  the  gamma  distribution  was  the  best,  four 
for  the  Welbull,  and  one  for  the  log-normal  (reference  table  6).  Generally,  Che 
gamma  distribution  was  the  best  for  the  wind  speed  distributions  during  all 
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seasons,  time  periods,  and  height  levels  with  5  exceptions.  The  best  distribu¬ 
tion  for  the  summer  night  period  for  height  level  3.0  Icm  was  the  log-normal.  The 
Welbull  distribution  was  best  for  the  summer  day  period  for  height  levels  0.3  and 
1.0  km,  annual  night  period  for  height  level  1.5  km,  and  annual  day  period  for 
height  level  0.5  km.  For  the  Mideast  the  oboerved  wind  speed  distributions  are 
more  "lumpy"  than  those  for  the  European  regions.  However,  the  analytical 
distributions  appear  to  adequately  represent  the  observed  distributions. 

6 .  CONCLUSIONS 

The  observed  wind  speed  distributions  are  adequately  represented  in  all  but  one 
case  by  the  Welbull  or  gamma  analytical  distributions  for  both  the  European  and 
Mideast  regions.  Generally,  the  Welbull  is  the  best  distribution  for  Europe  and 
the  gamma  the  best  distribution  for  the  Mideast.  One  might  use  the  analytical 
distributions  to  generate  wind  speeds  by  pseudo  random  numbers.  In  all  cases, 
generated  wind  speeds  in  excess  of  80  m/s  should  be  discarded.  Furthermore, 
occurrences  of  such  large  values  of  wind  speed  may  indicate  a  problem  with  the 
pseudo  random  number  generator. 
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Figure  I.  Regression  fitted  Weibuil  distribution  for  European  Lowlands  winter 
wind  speeds  for  height  level  0.5  )un. 
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Figure  3.  R*gr«8Ston  fitted  weibull  dietcibution  for  European  Lowlands  winter 
wind  speeds  for  heigfit  level  1.0  km. 
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Figur*  5.  Regression  ficced  Weibull  distribution  for  European  Lowlands  winter 
wind  speeds  for  neignc  level  l.S  )un. 
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Figur#  6.  R*grs>*i,on  gamma  disccxbution  for  European  Lowlands  wind 

speeds  for  height  ’evel  l.S  km. 
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Flgur*  8.  Ragrvsaion  fictad  gamma  di.scri.l3uti.on  for  Europaan  Lowlands  wind 
spaads  for  iiaigne  lavai  2.0  km. 
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Figur*  9.  RcgraaSLon  fittad  Waibuli  dutribution  for  European  Lowlands  wi.ncar 
wind  spaads  for  haighc  laval  2.S  )un. 
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Ragraasion  fictad  gamma  distri-bucion  for  Curopaan  Lowlands  wind 
spaads  for  haight  laval  3.S  km. 
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Figura  11.  R«gr«stion  {xccsd  w«ibull  diitrxbutxon  for  Curopacn  Lowland! 
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FigurB  12.  RcgrBition  fittad  gBrrun*  d].acribut.i.on  for  Europaan  Lowlands  wind 
spaada  for  haighc  laval  3.0  lun. 
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Figur*  14.  Regrasalon  flttad  gainina  distribution  for  Europasn  Lowlands  wind 
spasds  for  harght  laval  3.S  km. 
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Figur*  IS.  R€gr«8»ion  Jitttd  W«».buU  dutribution  for  European  Lowland* 
wuntac  wind  sptads  for  haighc  laval  4.0  lun. 
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Figura  16.  R«gr«aBi.on  fittad  waibull  distribution  for  Europaan  Lowianda  wind 
spaada  for  hai<5ht  laval  4.0  km. 
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Fi.t;uc«  20.  R«qr«««lon  fitted  w^ibull  distribution  for  European  Lowlands  wrnd 
spaads  for  haxght  iavvi  S.O  km. 


55 


09H 


Ul 


Figure  21 


R*gr*aaion 
wxnd  apaada 


EUROPEAN  HIGHLANDS 
WINTER  ALL  (LST) 

HEIGHT  LEVEL:  0.5  Km 

GAMMA  DISTRIBUTION  (REGRESSION) 

SHAPE:  2.04370 

SCALE:  2.29426 

NON-CALM  OCCURRENCES:  6759 

OBSERVED  DATA:  . 

FITTED  DISTRIBUTION:  - 

OBSERVED  -  FITTED:  . 


u» 

o 


ik 

Ul 


HIND  SPEED  (MPS) 

I  t.  I  -  i„ 


u> 

o 


WIND  SPEED  (MPS) 


4- 

Ul 


ad  gamma  di.8cri.bution  for  Eurcpaan  Hlghlanda  wi.n 
hai.ght  laval  O.S  km. 


09+  +  09H 


100  80  60  40  20  0  100  80  60  40 

PROS  (X)  CUM  PROB  (X) 


HIND  SPEED  (MPS) 


Figur*  22.  Regression  fitted  weibull  distribution  for  European  Highlands  wind 
speeds  for  height  level  O.S  km. 
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Flqurt  23.  R«gr«s*i.on  fitted  Weitjull  distribution  for  Europ««n  Highlands 
winter  wind  speeds  for  height  level  1.0  lun. 
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FLgur#  24.  R*gr*ssion  fitted  gamma  diatribution  for  Curopoan  Highlands  w^nd 
sp««ds  for  haight  Isval  i.o  km. 
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Figur«  25.  Ragraation  fictad  waibuil  diaccibucion  for  Europaan  Highlanda 
wvntar  wind  apaada  for  haignc  laval  1.5  km. 
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Figura  26.  R«gr*t»i.on  fitted  gaium*  ducr ibution  for  European  Highland*  wind 
spaada  for  haighc  Xaval  I.S  km. 
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Figuc#  27.  Mgr«*iion  fitted  diitribution  for  Europ*»n  Highl«ndi 

winc*r  wind  sp««ds  for  hdighc  2.0  kA. 
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Regrasaion  (Itcad  qamma  distribution  for  European  Highlands  wind 
speeds  for  hexght  level  2.0  km. 
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R«<}r«ation  fietad  waibull  diacribution  for  Europaan  Hiqblanda 
wintar  wind  apaada  for  haiqht  laval  2.S  kin. 
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Figur«  30.  R«gr«*«ion  fitted  q*mm«  distribution  for  EuropMn  Highlands  wind 
spssds  for  hsight  Isvsl  2.5  km. 
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Figure  32.  Regresalon  fitted  welbull  dletrtbutlon  for  European  Highland*  wind 
apeeda  for  height  level  3.0  kin. 
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fiqur#  33.  R«qr«««i.on  tittmd  Waibull  distribution  for  Europssn  Hiqhlsnds 
wintsr  wind  spssds  for  hsiqht  Isvsl  3.5  km. 
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EUROPEAN  HIGHLANDS 
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Figur*  34.  fitted  Wvibull  tbution  for  Eucop«»n  Highland*  wind 

tp««d*  for  height  l*v*l  3.S  iwi. 
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flgura  36.  Ragrataion  fltcad  WalDuil  distribution  for  Europaan  Highlands  wvnd 
spaads  for  hsight  laval  4.0  km. 
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Figure  37.  R«qr«»«i.on  fictad  W«i.bull  distribution  for  Curopcsn  HvghlAnda 
wintsr  wind  sp««da  for  haigtit  lavai  4.5  iun. 
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Figure  38.  Regreeeion  fitted  Wetbull  dietr ibution  for  European  Highland!  wj.nd 
speeda  for  height  level  4.S  )un. 
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Figur*  39.  R«qr«**ion  fitted  weibull  distribution  for  European  Highlands 
Winter  wind  speeds  for  height  level  5.0  kia. 
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Figur*  40.  R«gr«a«ion  fittad  Waibull  dlacribucion  for  Europaan  Highland!  wind 
spaada  for  haighr  laval  S.O  lun. 
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Figur*  41.  Regr«9Si.on  { icced  gamma  dxscr  ibuc  ion  for  Hi.d«aat  Oatarc  ni.ghttvma 
aummar  wind  apaada  for  haighc  laval  O.S  lun. 
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Figure  42.  Regrtssion  ticttd  WaibuXl  aiscribucton  foe  Hidaest  Oaaart  daytitna 
•urnmar  wind  spaada  foe  haeght  laval  O.S  km. 
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Figur®  43.  R«gr®*«i.on  <i.tt®d  ganvn*  distribution  for  Midlist  D®s«ct  nighttim® 
wind  sp««ds  for  height  l«v«i  O.S  km. 


78 


09  H 


PROB  (%)  CUM  PROB  (X) 


MIDEAST  DESERT 
ANNUAL  DAY  (LST) 

HEIGHT  LEVEL;  0 . 5  Km 

WEIBULL  distribution  (REGRESSION) 

SHAPE:  2.21271 

SCALE:  7.53769 

NON-CALM  OCCURRENCES:  3443 

OBSERVED  DATA:  . 

fitted  distribution:  - 

OBSERVED  -  FITTED: 


WIND  SPEED  (MPS) 


Figure  44.  Regression  fieeed  weibuil  distribution  for  Mideast  Desert  daytime 
wind  speeds  for  height  Level  O.S  kni. 
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Figur*  46.  R«gre*ii.on  fitted  Weibull  distribution  foe  Mideast  Desert  daytime 
summer  wind  speeds  for  berght  level  1.0  km. 
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Figur*  51.  R«gr«ifion  fitted  Weibuil  distribution  for  Midesst  Desert 
nighttime  wind  speeds  tor  height  level  1.5  kjn. 
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Figure  SJ.  Regression  fitted  gamm*  distribution  for  Mideast  Desert  nighttime 
summer  wind  speeds  for  height  level  2.0  km. 
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Figure  54.  F«gr*tBi.on  fitted  qamm*  di.str Lbut von  for  Mid«*«t  D«*«rt  daytim* 
wvnd  ipaada  for  h«vghc  l«v«l  2.0  km. 
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Figure  SS.  Regression  fitted  gairana  distribution  for  Hideast  Desert  nighttime 
wind  speeds  for  height  level  2.0  km. 
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Figur«  56.  Ragrassion  fitted  gamma  distribution  for  Hidsast  Oaaart  daytim* 
wind  apaada  for  haight  iavai  2.0  km. 
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figur*  S7.  KftgraasLon  exceed  gamma  distribution  for  Midaaat  Oasert  nighttime 
summer  wind  speeds  for  height  level  2.S  lun. 
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Pigur«  59.  Kcqrtation  {itcad  gamma  diacr;.buci.on  for  Mldaaat  Daaart  nrqhtt:.ma 
wind  apaada  for  haight  laval  2.5  km. 
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Figure  61.  degression  (icced  log~normel  distribution  for  Hidsast  Desert 
nighttime  summer  wind  speeds  for  height  level  3.0  pok. 


96 


100  00  60  40  20  0  100  80  60  40 

PROS  iX)  CUM  PROB  (X) 


-f- 


MIOEAST  DESERT 
SUMMER  DAY  (LST) 


Fiqur»  62.  fitted  dittrvbution  tor  Mi.d«««t  DoMrt  d*yci.m« 

tumrti«r  wind  •p««d«  for  htiijht  l«v«l  3.0  km. 


97 


MIDEAST  DESERT 
ANNUAL  NIGHT  (LST) 

HEIGHT  LEVEL:  3.0  Km 

GAMMA  DISTRIBUTION  (REGRESSION) 

SHAPE:  3.14780 

SCALE:  3.01649 

NON-CALM  OCCURRENCES:  8949 

OBSERVED  DATA:  . 

FITTED  DISTRIBUTION:  - 

OBSERVED  -  FITTED:  . 


4- 


- - 1 - -  I 

UJ 

o 

WIND  SPEED  (MPS) 

_ » _ i _ _ _ 


Ik. 

cn 


T 


mideast  desert 

ANNUAL  DAY  (LST) 

HEIGHT  LEVEL:  3.0  Km 

GAMMA  DISTRIBUTION  (REGRESSION) 

SHAPE:  3.12657 

SCALE:  2.91417 

NON-CALM  OCCURRENCES:  11658 

OBSERVED  DATA:  . 

FITTED  DISTRIBUTION:  - 

OBSERVED  -  FITTED:  . 


“h 

CJI 


- - - - h 

o  ^ 

WIND  SPEED  (MPS) 


01 

o 


R«<jr«ttion  fitted  qamma  dltcr ibutLon  for  Midaaat  Oaaart  daytima 
wind  apaads  for  haight  laval  3.0  km. 
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Fi.gur«  67.  R«gr««*ion  fitted  gamma  <aitcri.buti.on  for  Kldaatt  Oatart  niqntti,ma 
wind  tpaada  for  havght  lavai  3.5  km. 
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figuce  63.  Regr«S8ion  fictcd  tjamma  dLateibuti.&!\  for  MidtiSt  Dasact  daytiirts 
wind  spaads  for  heighc  laval  3.3  !un. 
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figure  69.  Regreaaion  fitted  ganwa  distribution  foe  Mideaet  Desert  nightt:.T,e 
summer  wind  speeds  for  height  ievel  4.0  Ion. 
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Figur*  70.  Regr»f«ion  fitted  gamma  diatribution  for  Mideaat  Desert  daytime 
summer  wind  speeds  for  height  Level  4.0  km. 
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Figure  12.  Regression  fitted  gemms  distribution  for  Mideest  Desert  dayt4.me 
Mind  speeds  for  height  level  4.0  kin. 
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Fi9ura  73.  Regression  fitted  gamnte  dietridution  for  Hidaaat  Desert  nighttin^.e 
Summer  wind  speeds  for  height  level  4.5  km. 
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Figur#  74.  Regression  fitted  gamme  distribution  for  Mideast  Desert  daytime 
SurTwer  wind  speeds  for  height  level  4.S  kin. 
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Figure  75.  Regression  fitted  gemma  distribution  for  Mideast  Desert  nighttime 
wind  speeds  for  height  level  4.S  km. 
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Figure  77.  Regreeeiort  fitted  gemma  distribution  for  Mideaat  Detect  nighttime 
summer  wind  speeds  for  height  level  5.0  )«n. 
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Pigure  73.  Regression  ficted  gamma  distribution  for  Mideast  Desert  daytime 
summer  wind  speeds  for  height  level  S.O  km. 
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Figur#  79.  Ragseasion  fitted  gamma  distribution  for  Midaast  Dasart  nighttime 
wind  speeds  for  haight  level  5,0  km- 
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Figure  80.  Regreaaion  fitted  gemma  dietribution  foe  Hideaet  Deeart  daytime 
wind  speeds  for  height  Level  S.O  km. 
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